Non-Rayleigh limit of the Lorenz-Mie solution and suppression of the scattering 

by the spheres of negative refractive index 
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We study light scattering by spheres made of negative refractive index materials. We demonstrate 
that scattering efficiency of the light does not always obey the well-known Rayleigh's law Qsca ~ 1/A* 
for small values of size parameter, but could be wavelength independent Qsca ~ const, or even 
increase with the wavelength as Qsca ~ A^. In addition, it is found that spheres with radius 
comparable with the input wavelength may exhibit suppression of the scattering in any directions, 
particularly in forward one. It can be considered as analogues of Brewster's angles for spheres. 

PACS numbers: 42.25.Bs, 42.68.Mj, 92.60.Ta 



I. INTRODUCTION 



Light scattering by small particles is one of the fun- 
damental problems, which found a broad application in 
various fields such as meteorology, biology and medicine. 
The first results on scattering by small particles were de- 
rived by Lord Rayleigh 0, [3, S], and analyzed by J.C. 
Maxwell-Garnett 0, 1^]. It is widely excepted that the 
complete solution for the homogeneous spheres of arbi- 
trary size and optical constants was found by Mie more 
than a century ago @. As a matter of fact, the first ex- 
act solution was obtained by Lorenz in 1890 [Tj], based 
on his own version of the electromagnetism. But, this 
paper was published in Danish, and did not receive a 
proper acknowledgement by the scientific community Q . 
Thus, it will be reasonable from a historical point of 
view to call it Lorenz-Mie solution. This exact solu- 
tion was carefully analyzed, and successfully applied to 
describe scattering by particles of various shapes. In 
the limit of small spheres it reproduces the results ob- 
tained by Lord Rayleigh. However, recent studies of res- 
onant scattering by small particles with weak dissipation 
rates have revealed novel and unexpected features, 

namely giant optical resonances with an inverse hierar- 
chy (the quadrupole resonance is much stronger than the 
dipole one, eta) [Hi , a complicated near- field structure 
with vortices |12l | , unusual frequency and size dependen- 
cies [isj . which allow to name such a scattering anoma- 
lous. 

Nowadays technological success in fabrication of nanos- 
tructured materials has led to an emergence of the rapidly 
growing field of nanophotonics. The primary interest of 
this field is to achieve an enhancement of light-matter in- 
teraction at the nanoscale. The potential for this is based 
on the possibility to excite localized surface plasmons in 
metallic nanoparticles. The strong light confinement in 
ultra-small volumes may result in various nonlinear ef- 
fects. The Lorenz-Mie solution is able to shed a light on 
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the origin of such enhancement by metal nanoparticles, 
and to provide better understanding of the phenomenon. 

The beginning of XXI century was marked by open- 
ing a new chapter of electrodyanical theory, when it was 
generally accepted by scientific community that material 
with negative refractive index can exist, and can even be 
fabricated. A lot of attention was attracted to such kind 
of materials with negative index of refraction since they 
possess unusual optical properties such as inverted Snell's 
law [13], backward waves propagation, possibility to de- 
sign perfect lens [31, and even cloak of invisibility [3. 
The origin of such peculiar effects lies in the interplay of 
dielectric and magnetic responses. Almost every corner 
of electrodynamic theory was revised from the perspec- 
tive of materials with negative refractive index with aim 
to find novel phenomena, which were overlooked in the 
past due to restricted parameters range. 

In this paper we study light scattering by spheres made 
of materials with negative refractive index. For that 
purpose we are using exact Lorenz-Mie solution, which 
is valid for any values of permettivities and permi- 
abilities /ig, including negative ones, and various ratios 
K = 2-Ka/X of radius of the sphere a to the wavelength 
A. This solution predicts that spheres of the small radius 
a — > made of conventional materials will scatter light 
as Qsca ~ l/-'^^, which is in agreement with the Rayleigh 
scattering by small particles. It means that the blue com- 
ponent of the light is scattered 10 times stronger than 
the red one. We evidence this effect everyday by looking 
at the blue sky. Another prediction of Lorenz-Mie solu- 
tion for conventional materials is that scattering by large 
spheres is predominantly in forward direction, which is 
knows as Mie scattering. According to our analysis of 
the Lorenz-Mie solution it turns out that both effects 
(Rayleigh and Mie scatterings) can be altered by nega- 
tive refractive index materials. In particular, the light 
scattering by small spheres can be wavelength indepen- 
dent Qsca ^ const, or even diverge Qsca A^ . In addition 
to this, the scattering can be significantly suppressed in 
any direction, including the forward one. 

The paper is organized as follows. In Section II the 
Lorenz-Mie solution is described. The scattering by small 
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sphere with negative refractive index is studied in Sec. 
Ill, where different regimes are identified. In Sec. IV 
the angular dependence of scattered light is considered. 
Section V concludes the paper. 



II. LORENZ-MIE SOLUTION 

We consider the problem of light scattering by a spher- 
ical particle of radius a, made of an material with arbi- 
trary permettivity es and permiability /Zs , embedded into 
surrounding media with given optical constants Eq and 
/io. For that purpose we will utilize the exact Lorenz-Mie 
solution ■ According to this solution the scattering 

efficiency of light can be written as follows 



r.2 



(2n + l)(|a„|2 + |6„|2) 



(1) 



where n = lnajX is the size parameter, a„ and 6„ are 
scattering coefficients. The latter ones can be expressed 
in terms of Riccati-Bessel functions of the first and third 
kind 



(2) 



where m — \/e^Ji7/£ol^ is the ratio of refractive indices, 
and r] = mjjLo/ Us = \J Esy^ol^o^is is the ratio of elec- 
tromagnetic impedances of the sphere and surrounding 
media. Due to the symmetry of the scattering coeffi- 
cients [is'] the expression for 6„ coefficients can be ob- 
tained as follows 



bn{r],m,K) = a„(l/?7,TO, k) . 



(3) 



In the case of surrounding media with equal optical con- 
stants Co = fJ'o (air, for example), there are other sym- 
metries of the scattering coefficients, which will be used 
later. When the permettivity and pemiability of the 
sphere are the same eg = ^is, the parameter -q reduces 
to unity. Due to relation ^ the scattering coefficients 
become also identical a„ — bn- The substitution — > jig 
and ^s — * Es results in 77 — > l/?7, which again due to rela- 
tion ^ will simply swap scattering coefficients a„ &„ 
and hn — > a„. Both symmetries leave the scattering ef- 
ficiency ([1]) unchanged. Because of these symmetries, it 
will be enough to consider materials, whose optical prop- 
erties satisfy the relation |es| < \ns\- 

In the case of nonabsorbing sphere, /m(es) — 
Im(fis) — 0, the scattering efficiency is equal to the ex- 
tinction efficiency Qsca = Qext, and, thus, describes the 
overall attenuation of the light as it propagates through 
the sphere. 



The angular dependence of the scattered light can be 
studied by using scattering amplitudes 

00 2n + 1 
= , -, j Qn7r„ + b„T„) , (4) 
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of perpendicular Sj_ and parallel S\ | polarized light to the 
scattered plane, determined by the incident and scattered 
directions. The functions 7r„ and r„ can be defined via 
Legendre polynomials 
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From expressions Q we can define transmission T (for- 
ward, e = 0°) and reflection R (backward, 6 = 180°) 
amplitudes 

T = 5i,||(0°) = -^(2n+l)(a„ + 6„), (6) 

n=l 
1 °° 

R = ±5^,11(180°) = -^(-I)"(2n+l)(a„-6„) . 



n=l 



These two directions are special ones, since the scat- 
tering does not depend on the input polarization. 



III. (NON) RAYLEIGH LIMIT 

For spheres with the radius much smaller than wave- 
length of the light K ^ I the general solution ((T|) con- 
verges very fast, and it is possible to retain only the first 
terms ai and 61, since all other are negligible |a„| w 
6„| ^ for n > 1. The Riccati-Bessel functions for 
n = 1 have a very simple form [l3| 



V'i(p) = ^-cos(p) , ei(p) 



exp(*p) (- + !).(?) 
P 



To make our analysis applicable to any values of optical 
constants and /ig we will use relations ([7]) for functions 
with argument nin, and their expansions up to the second 
order [IS] 
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for functions with argument k in expressions ([2|). Thus, 
after some algebra the scattering coefficient ai can be 
written in the following form 
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Al 2k^{[1 + 2?7m][TOa;cos(mK) — sin(TOK)] + TO^K^ sin(mK)} 

A2 rim[6i — 3zk^ + 4^^] [rriK cos(mK) — sin(mK)] + [2^^ — 6i — 3iK^] [rriK cos(mK) — (1 — rn^n^) sin(mK)] 



The expression for bi scattering coefRcient can be ob- 
tained by using the symmetry relation ([3]). 

To make a proper analysis of the scattering coefficient 
ai in the limit of the small size parameter k 1, we ex- 
pand the numerator Al and denominator A2 separately 
up to the same order 



Al 



A2 = -2i{'qm + 2)m^ 



(10) 
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10 -|- rim{5 + m )]m k' 



^{rim-l)m^K^ + 0{k'^) . 

In general only leading orders are kept in Al and A2, 
which allow one to simplify the scattering coefficients ai 
and bi to the form [l^ 
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The substitution of these relations into the expression of 
the scattering efficiency ([T]) will result in 



Qs 





e-s - Eo 


2 

+ 




- Mo 


2" 




Us + 2eo 




Ms - 


h2/io 





(12) 



which is in agreement with Rayleigh's approach, accord- 
ing to which small spheres effectively scatter light as 
Qsca ^ l/A**. This is a well-known result which is appli- 
cable to spheres made of almost any materials, except for 
Es = — 2eo or /is = — 2/io. In this case one of the scatter- 
ing coefficient diverges ([TT|) , and so does the scattering 
efficiency (IT^ . Usually, this divergence is associated with 
the excitation of a localized surface modes of spheres of 
the larger radius k > 1 17]. But, this interpretation is 
not applicable in the limit of small spheres k <^ 1, and 
should be revised. 

Indeed, the relations (fTT|) are only expansions of the 
scattering coefficients oi and hi And if an expansion 
produces some divergences we should look back into the 
original expression and take a special care for these par- 
ticular parameters. In the expressions (|10p the first term 
of the denominator A2 vanishes exactly when Cs = —2eo- 
It means that we should take the second term to obtain 
the correct expression for the scattering coefficient oi 
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FIG. 1: (Color online) Loglog plot of the scattering efficiency 
dependence vs size parameter k, calculated by using the ex- 
act Lorenz-Mie solutions for the sphere embedded into air 
(to = Mo = !)■ We chose three sets of parameters to illus- 
trate different types of the asymptotic behavour : i) = 2 
and /is = 1 (dotted line); ii) es = —2 and /is = —3 (dashed 
line), iii) es = —2 and /is — —5 (solid line). The calcu- 
lated results are in very good agreement with the analytical 
ones (|12I14I16|I . Numerical errors in calculations of the Bessel 
functions become essential for the very small size parameter 

K < 10"^ 



while the leading order of bi scattering coefficient remains 
unchanged. The scattering efficiency in this case becomes 
wavelength independent 
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One may notice that if in addition to the condition eg = 
— 2eo the permeability of the sphere satisfies the relation 
fig — ~5fio both the scattering coefficient oi p3p and 
efficiency Qsca p4|) will diverge. It means that one should 
revise the expression for the denominator A2 again. For 
these particular parameters first two terms vanish in A2 
(fTU)) . and the third term is identical to the nominator Al. 
It results in the following scattering coefficients 



ai « 1 , bi 



Ai 



(15) 



The scattering efficiency in this case can be written as 
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(16) 
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(a) Electi'ic field distribution |Ep (b) Magnetic field distribution |M| 




FIG. 2: (Color online) The distribution of the electric (a) and 
magnetic (b) fields inside the sphere of radius a = 10~^. 



which is proportional to the square of the wavelength 
Qsca ^ A^. We would hke to note here, that the expres- 
sion similar to eq. ([9]) was obtained in Ref. [l^l without 
subsequent analysis. We found it incorrect since it does 
not reproduce the Rayleigh scattering regime, although 
it still allows to predict correctly other limits. 

The scattering efficiencies calculated by using the ex- 
act Lorenz-Mie solution are shown in Fig[TJ The asymp- 
totic behaviour of all curves are in a very good agree- 
ment with analytical derivations (|12ll4ll6p . except the 
very small values of the size parameter x < 10~^ because 
of the difficulties of numerical calculations of the Bessel 
functions \1T\ . 

Naturally, the question arises about the origin of 
changes of the scattering efficiency asymptotic behaviour. 
As it was mentioned above, the argument of localized 
surface mode excitation is not applicable here, because 
of the extremely small radius of the sphere compared to 
the wavelength. To answer this question, we calculate 
the distribution of the electric and magnetic filcds inside 
the sphere. It turns out that the distribution of the elec- 
tric field is the same for all parameters, and resembles 
the electric dipole excitation [see Figl^Ja)], while the in- 
tensity scales with actual values of the optical constants. 
But what is changing inside the sphere is the distribu- 
tion of the magnetic field. When eg ^ — 2eo the magnetic 
field inside the sphere is almost homogeneous and prac- 
tically zero. The drastic change of the magnetic field 
distribution takes place when the condition = — 2eo 
is satisfied. It is no longer homogeneous, and resembles 
magnetic dipole excitation [see Figl^l^b)]. In general, un- 
der this condition, the magnetic field is still very weak, 
but becomes order of unity when the additional condition 
Us = —5^J'o is satisfied. Thus, it allows us to conclude 
that the excitation of the magnetic dipole is responsi- 
ble for the variation of the asymptotic behaviour of the 
scattering efficiency. 



IV. SUPPRESSION OF THE SCATTERING 
EFFICIENCY 

The analysis of the angular dependence of the scat- 
tered light ^ suggests that the scattering in the back- 
ward direction ([6]) can be totally suppressed i? = when 
the scattering coefficients are the same a„ = &„ for all n. 
This condition is realized when the optical properties of 
the sphere are identical — (UlEl), due to the sym- 
metry of the electric and magnetic fields. On contraray, 
the scattering in the forward direction ^ can be sup- 
pressed r = 0, when the scattering coefficients are oppo- 
site a„ = —bn for all n. Note here, that in both cases the 
strength of the overall scattering efficiency is the same 
^ , indicating the fact that the sphere is equally excited 
under both conditions. Obviously, they can not be ex- 
actly satisfied simultaneously, otherwise a„ = 6„ = 0. 

As it was shown by Kerker [T^ suppression in the 
forward direction can be achieved for small magnetic 
spheres. Indeed, according to (|lip for magnetic spheres 
embeded into air the condition ai = — &i is fulfilled for 



In Ref. [2l| it was demonstrated that there is one ex- 
ception of this formula, namely =^ fi^ = —2, which 
does not produce zero scattering in the forward direc- 
tion. This is a very interesting case, because at the first 
sight it appears that both conditions for zero scatter- 
ing in forward ai — — 6i [because of p7)) ] and back- 
ward ai — bi [because of es = ^j,s] directions are sat- 
isfied, which should require that both scattering coeffi- 
cient should vanish oi = bi = 0, s ugg esting that sphere 
is transparent, which is not true [2l|. But the answer 
to this paradox is quite simple. Our analysis in the pre- 
vious section of the scattering coefficients for the very 
small spheres suggests that expansions (ITT|) are no longer 
applicable because of divergence at Cs — —2, and the 
proper expansions should be used instead, accord- 
ing to which ai ^ —bi. As a result, the relation ^T7\ is 
not applicable in this case, and, therefore, the scattering 
in the forward direction is nonzero. 

It is known that for spheres of large radius the scat- 
tering occurs predominantly in the forward direction 
T > 'S'j,jj(6' ^ 0°) 17]. But, recently we has demon- 
strated [2J| that weakly absorbing particles can exhibit 
suppression of the scattering in the forward direction, 
which can be explained in terms of the Fano resonance. 
This analysis is quite general, and can be extended to 
magnetic spheres as well. In Figure [Sja) we plot the de- 
pendence of the forward scattering T versus size param- 
eter K for optical constants es = /i^ = — 5. It indicates 
that near the size parameter value k « 0.827 there is 
a resonant suppression of the forward scattering T — > 0. 
According to the optical theorem |17j | , which can be writ- 
ten in the absorptionless case as 

Qext = Qsca = —ReiT) , (18) 
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(b) K = 0.827 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

size parameter, k 

FIG. 3: (Color online) The scattering by the sphere with opti- 
cal constants es — fis = —5: (a) the scattering efficiency Qscn 
(dashed line) and the forward scattering |Tp (solid line) vs 
size parameter k, (b) angular dependence of the light scatter- 
ing in the vicinity of the suppression of the forward scattering 
at K « 0.827. In this case the scattering is predominantly in 
the transverse directions. 



the suppression of the scattering in the forward direction 
will result in the overall suppression of the extinction of 
light, making the sphere effectively transparent. Indeed, 
Fig. [2{a) clearly demonstrates that there is suppression 
of the scattering efficiency in the vicinity of the suppres- 
sion in the forward direction. It is interesting to look 
at the angular dependence of the scattered light, since 
in this particular case the condition for zero backward 
scattering for any value of the size parameter is satisfied 
Es = A^s • The resulting plot is showing in Fig. ^h) . It re- 
sembles figure-8 profile. It means, that all light is mostly 
scattered in transverse directions only. This is quite un- 
usual behaviour, and it is opposite to the Rayleigh angu- 
lar scattering, which is 90° rotated, and resembles figure- 
oo. 

This analysis can be extended to the scattering in any 
direction 0° < 6* < 180°. In general, each polarization 
S'xji will be scattered differently In case of magnetic 
particles it is possible to find the conditions when one of 
the polarizations will be strongly suppressed, resulting in 
totally polarized scattered light. These conditions can be 
considered as generalized Brewster's angles of the sphere. 
In case of the sphere of the small radius, such conditions 
can be found analytically from (|4|) and ([6]) with tti = 1 
and Ti = cos 9 



cose* 



1 



II 



according to which S^{6*^) = and S'|H,i/||y 
complete analysis of the angular depended scattering by 
small magnetic spheres can be found in Ref . [T9l| . 

For spheres with the radius larger or comparable with 
the wavelength we can perform numerical analysis in any 
direction 9* , similar to the forward scattering in Fig.[3ja). 
The dependences of the scattering amplitudes ^ ver- 
sus size parameter k for optical constants = — 1 and 
/is = —5 in the direction 9* = 45° are presented in 



(19) 



0. The 



(a) 


1 ^^S|(45"||! 






|S^(45-)|! 



(b) 









,815 1.S2 



size parameter, K 



size parameter, K 



FIG. 4: (Color online) The scattering by the sphere with op- 
tical constants = — 1 and /is = —5 in 6* = 45° direction 
vs size parameter k: (a) scattered perpendicular S± (dashed 
line) and parallel S\\ (solid line) polarized light; (b) total po- 
larization P according to ea. (|20|) . 



Fig. Ua) . This particular example illustrates that when 
one of the scattering components iS'_l is resonantly sup- 
pressed another one 5*1 1 is resonantly enhanced in the 
vicinity of k « 1.8115. In addition, there is resonant 
suppression of S\\ at k « 1.81. All of these results in a 
sharp variation of the polarization 



15, 



l^lll^ 



(20) 



The scattered light is almost perpendicular polarized 
\S±\'^ ~ 1 at K « 1.81, and becomes totally parallel po- 
larized IS'iiP = 1 at /« 1.8115 [see Fig.|4i;b)]. In other 
words, by changing slightly the input wavelength, the 
scattered light in 9* = 45° will swap it polarization. One 
might notice an asymmetric profile of resonant shapes, 
which, actually, can be understood in terms of the Fano 
resonance [221. 



V. CONCLUSIONS 

We have studied light scattering by spheres made of 
materials with negative refractive index. By using the 
exact Lorenz-Mie solution we have demonstrated that 
in addition to the well-known wavelength dependence of 
the scattering efficiency Qsca ~ 1/^"'', the negative in- 
dex materials may inverse this dependence Qsca ~ A^, 
or even lead to the wavelength independent behaviour 
Qsca ~ const. The standard wavelength dependence was 
obtained by Lord Rayleigh in the assumption that small 
sphere can be represented by an electric dipole. Our 
analysis suggests that for certain values of the optical pa- 
rameters Cs and /is of the sphere a magnetic dipole may 
become excited as well. Its excitation results in a qualita- 
tively different wavelength dependence of the scattering 
efficiency. Such regions have been found analytically, and 
confirmed numerically by using the Lorenz-Mie solution. 
We have also analyzed angular dependence of the scat- 
tered light for spheres with the radius comparable with 
the wavelength. It was found that both forward and 
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backward scattering can be effectively suppressed lead- 
ing the scattering predominantly in the transverse direc- 
tions. Recent experimental observations [13, [13, [IE, [111 
of the suppression of the forward scattering by magnetic 
spheres are in agreement with our results. Moreover, we 
have demonstrated that the suppression of the scattering 
can be achieved in any given direction for various input 
polarization, providing with totally polarized scattered 
light. In analogous with the Brewster's angle at the in- 
terface of two media, such angles can be called Brewster's 
angles of the sphere. 

It is worthwhile noting that the prefactor (cs — 
em)/(es + 2em) in expression of the scattering coefficient 
ai (fTTj) is exactly the same as Clausius-Mossotti factor in 
Lorenz-Lorentz formula, describing the relation between 
the refractive indices and polarizability. Our analysis 
suggests that in the case of negative refractive index ma- 



terials, and in the limiting cases considered above, the 
Lorenz-Lorentz formula might be also changed according 
to eqs. (|13|15p . which will require a separate study. 

Thus, if the success in fabrication of metamaterials 
with negative refractive indices will be achieved we will 
be able, after all, to alter the familiar colour of the sky 
to be red, or even white. 
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